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The quantum states and thermodynamical properties of the Te vacancy in CdTe are addressed 
by first principles calculations, including the supercell size and quasiparticle corrections. It is shown 
that the 64-atoms supercell calculation is not suitable to model the band structure of the isolated Te 
vacancy. This problem can be solved with a larger 216-atoms supercell, where the band structure 
of the defect seems to be a perturbation of that of the perfect crystal. It is interesting to note 
that the Te-vacancy formation energy calculated with both supercell sizes are close in energy, which 
is attributed to error cancelation. We also show that the interplay between supercell size effects 
and the band gap underestimation of the generalized gradient approximation strongly influences the 
predicted symmetry of some charge states. 
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I. INTRODUCTION 

Cadmium telluride (CdTe) is a semiconductor crys¬ 
tal employed as photovoltaic material which is growing 
rapidly in acceptance, becoming a cheap alternative to 
the silicon solar cells. CdTe has advantages for photo¬ 
voltaic applications, such as a direct band gap, allowing 
to absorb solar energy in thin films. Its band gap of 
1.5 eV at room temperature is in the range for maximum 
conversion efficiency of a single absorber solar cell.i CdTe 
is also widely used for radiation detectors;^ However, the 
nature of the defects in CdTe thin films has been sub¬ 
ject of controversy for long time, and its identification 
frequently relies on thermodynamics arguments and ab 
initio calculation of the defects and impurities formation 
energies. 

Electron paramagnetic resonance (EPR) spectra un¬ 
veiled a paramagnetic center that was initially assigned 
to the tellurium vacancy Initial evidences suggested 
that its charge states are neutral and singly positive, with 
the transition level e(-|-/0) = 0.2 eV^® However, this as¬ 
signment has been rejected based on both density func¬ 
tional theory (DFT) calculations of total energies^ and 
hyperfine interaction termsj^ Later on, this EPR signal 
was reassigned to the Tecd antisite<i No direct evidence 
of Vxe bas been obtained afterwards from experiments, 
but from theoretical calculations this defect is one of the 
most important intrinsic donors 

The formation energy of Vxe can be obtained com¬ 
puting the difference between the total energies of a 
crystal supercell containing n CdTe units, i?(Cd„Te„), 
and the energy of the same supercell with a missing Te, 
i?(Cd„Te„_i), supplemented with the Te chemical po¬ 
tential £'(Te) -I- A/TTe 

AiT/(VTe^) = E(Cd„Te„_i) + E(Te) - E(Cd„Te„) 

+A/rTe + q[Ev + Ep] -\- A.Esize + ^Eq p . (1) 
In Eq.m i?(Te) is the energy of Te in a reference state, in 


this case it is the energy per atom of bulk Te. A/xxe is the 
chemical potential relative to the reference state, and it 
depends on the thermodynamic equilibrium conditions. 
In Te-rich conditions A/xxe = 0, while in Te-poor con¬ 
ditions A/ixe equals the CdTe formation energyFor a 
charged vacancy, q is the number of electrons donated to 
the environment, and Ey -\- Ep is the electron chemical 
potential expressed as the sum of the valence band max¬ 
imum Ey and the Fermi level Ep. AEgize and AE^.p. 
are size and quasiparticle corrections that are explained 
in appendices El and [B1 respectively. 


Most calculations of CdTe point defects have been in¬ 
vestigated using the 64-atom supercell (SC64) of the per¬ 
fect crystali^i^^— This supercell is 2 x 2 x 2 times the 
8-atoms cubic unit cell. In this article, we show that 
the SC64 supercell is too small to obtain reliable re¬ 
sults for the Vxe band structure, while thermodynamic 
calculations are hampered by uncontrollable size effects 
that are incompatible with the notion of diluted defects. 
We observe that changes in the conduction band are so 
dramatic that destroy the reference to apply quasiparti¬ 
cle gap corrections to the formation energy of negatively 
charged states. We show that this problem is solved us¬ 
ing a 216-atom supercell of the perfect crystal (SC216), 
which is 3 X 3 X 3 times the 8-atom cubic unit cell. In 
addition, the band filling induced by the use of a k-point 
grid causes a spurious distortion in the double negatively 
charged state 2—. We show that this size effect can be 
avoided by analyzing the band diagram and applying the 
band-filling correction. According to our results, a large 
supercell with only one k-point is enough to avoid this 
problem. For CdTe, this large supercell can be 4 x 4 x 4 
times the cubic unit cell, that contains 512 atoms for the 
perfect crystal (SC512). 
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II. COMPUTATIONAL DETAILS 

Our DFT calculations were performed using a plane- 
wave projector augmented wav o^^d^ scheme, as imple¬ 
mented in the Vienna Ab Initio Simulation Package 
(VASP)ii^ We have used mostly the generalized gradi¬ 
ent approximation (GGA) for the exchange-correlation 
functional as proposed by Perdew, Burke, and Ernzer- 
hof (PBE) In certain cases, we have also used the hy¬ 
brid functional proposed by Heyd, Scuseria and Ernzer- 
hof (HSE06)ii^ This functional generally allows one to 
obtain better band gap energies and structural proper¬ 
ties than those obtained with PBE, but at the cost of a 
great increment in computer time. 

We have used a plane-wave cutoff energy of 460 eV for 
the wavefunctions. To sample the reciprocal space, we 
have used T-centered 3x3x3 and 2x2x2 k-point grids 
for the SC64 and SG216 supercells, respectively. This 
setup allows convergences of 0.5 meV/atom and 0.4 kbar 
in total energy and pressure, respectively, as tested for 
variable cell relaxations of the GdTe zincblende unit cell 
with 8x8x8 k-point grid. In fact, a reduced cutoff of 
285 eV allows a convergence within 10 meV/atom in total 
energy and 2 kbar in pressure, but the properties derived 
from relative energies, e.g., lattice constants and bulk 
moduli are well converged. The GdTe lattice constant 
obtained with variable cell relaxation with a cutoff of 285 
eV differs by 0.4 % from the converged value, but when 
the energy is fitted with the Birch-Murnaghan equation 
of state, the difference falls to 0.005 %. The energy dif¬ 
ference between the zincblende and wurtzite structures 
is 3.5 meV/atom considering both cutoffs. The nxnxn 
supercell size is n times the converged theoretical lattice 
constant of the perfect crystal (6.6212 A). We have used 
a cutoff of 285 eV to relax the point defect structures 
with fixed supercell size, and a cutoff of 460 eV to com¬ 
pute the total energies. The highest precision is needed 
because the energy error is amplified by the number of 
atoms in the supercell, e.g., for the neutral Cd vacancy 
in SC64 supercell, the formation energy computed with 
both cutoffs differs by 0.03 eV. We have verified that for 
higher cutoff energy, the forces on the atoms in the re¬ 
laxed geometry keep below 0.01 eV/A. 

In previous works, the model SG64 has generally been 
used with 2x2x2 T-centered k-point gridsi^^— This 
sampling is equivalent to using only the T-point for the 
SG512 supercell. Eor , using the 2x2x2 and 3x3x3 
grids, the difference in formation energy is 0.030 eV for 
model SC64. For model SG512 with F-point sampling, 
the error in total energy is scaled to 0.24 eV. 

The calculated GdTe heat of formation is found to be 
of 0.93 eV, in good agreement with the available exper¬ 
imental value of 0.96 eVi^ The energy of bulk Gd has 
been obtained by relaxing the 2-atom unit cell of the 
experimental structure IGSD ^52793 available from the 
MedeA® software environmentThe Brillouin zone was 
sampled by a 22 x 22 x 10 k-point grid. The energy of 
bulk Te has been obtained relaxing the crystallographic 


3-atoms unit cell IGSD ^23060. The Brillouin zone was 
sampled by a 8 x 8 x 5 k-point grid and tested with a 
15 X 15 X 10 grid. 

The formation energies have been corrected for the ef¬ 
fects of supercell size and quasiparticle corrections to the 
band edges, as described in the Appendices. 


III. RESULTS AND DISCUSSIONS 

A. The Te vacancy in a 64-atom supercell: 

Structural and electronic properties 

The SG64 supercell is the standard framework nowa¬ 
days to model defect in diamond- or zincblende-type 
semiconductors. It is a popular choice because of its 
cubic shape, that maximizes the distance between the 
defect and its periodic images, and allows easy charge 
corrections. In addition, it is affordable for high level 
methods beyond GGA of DFT. The next cubic super¬ 
cell, containing 216 atoms, is still extremely costly for 
high level methods. Next we analyze the band struc¬ 
ture of the defective supercells, and show the problems 
to apply band-filling and quasiparticle corrections with 
the SC64 model. 

Figure [1] shows the band structure calculations of the 
perfect GdTe and Vxe in the SC64 supercell. The atomic 
positions were relaxed keeping the symmetry, while 
distorted configurations will be considered separately us¬ 
ing larger supercells. The vacancy in the SC64 supercell 
is strictly a periodic model of Vxe in a concentration of 
one every 64 atoms or 10^° cm“^. In Figure[T]we see that 
for the charge states 0,1—, and 2—, the conduction band 
(CB) is dramatically modified with a significative incre¬ 
ment of the band gap. This global change produced by 
the vacancy is incompatible with the notion that an im¬ 
purity should introduce small perturbations on the crys¬ 
tal band structure. This example shows that band dia¬ 
grams, which are in principle not adequate to study ape¬ 
riodic systems, are still useful to uncover flaws in defect 
models. From the practical side of thermodynamic cal¬ 
culations, the drastic alteration of the conduction band 
makes impossible to apply gap corrections, resulting in a 
false stability of the negatively charged states in n-type 
doping conditions. The only solution for the exaggerated 
influence of Vxe is to explore configurations with a larger 
supercell. 


B. The Te vacancy in the 216-atom supercell: 

Structural and electronic properties 

Figured] shows the band diagram of the perfect crystal 
and the Vxe defect in different charge states computed 
with the 3x3x3 supercell of the conventional 8-atoms 
unit cell (SC216). The Tj, symmetry of the atomic ar¬ 
rangement around the vacancy has been conserved in the 
relaxation as for SC64. Distorted configurations with 
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FIG. 1. Band structure calculations for the Te vacancy in CdTe performed with a 64-atom supercell, (a) Perfect crystal, (b) 
to (f) Te vacancy in the charge states 2 + to 2—. The dashed line indicates the highest occupied level. 




FIG. 2. Band structure calculations for the Te vacancy in CdTe performed with a 216-atom supercell, (a) Perfect crystal, (b) 
to (f) Te vacancy in the charge states 2 + to 2—. The dashed line indicates the highest occupied level. 


lower energy will be studied separately below. For all 
the charge states considered, the effects of the vacancy 
are small band shifts and splits. In addition, in Fig. [2] 
we can recognize an empty-defect band over 3.2 eV for 
g = 2—,1 —,0. As we will see later, the defect band 
becomes active, driving the so called Jahn-Teller distor¬ 
tions. The origin of the defect band can be understood 
as follows: Each Te-Cd covalent bond is populated by 0.5 


Cd electron and 1.5 Te electron. When Vxe is formed, 
four sp^ bonding orbitals and six electrons are removed. 
This process leaves two electrons from the four surround¬ 
ing Cd 5s orbitals. These orbitals combine to form four 
defect bands: one mixed with the host valence band (VB) 
that gets occupied by the electrons left by the Cd broken 
bonds, and three empty isolated bands (~ 3.1 — 3.5 eV). 
Of these three bands, the lowest two are degenerate. In 
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FIG. 3. (Color online) Stable geometries and band diagrams 
for the nentral Te vacancy withont (a) and with (b) a C 2 v 
distortion, nsing the SC216 snpercell. The solid black (green) 
lines are the bands of the defective (perfect) supercell, and 
the dashed lines represent the highest occupied level. Pink 
and cyan balls represent Te and Cd, respectively. 


states 1— and 2—, the defect levels keep empty and the 
extra electrons occupy perturbed host states (PHS)f22 
that in the supercell calculation appear similar to the 
CB. 

The neutral state shows a similar band structure if the 
atomic positions keep the Td symmetry. However, the 
total energy and the band diagram are modified due to 
a symmetry break, inducing a C 2 v distortion that low¬ 
ers the total energy by 0.3 eV. The distortion leads to 
formation of a Cd 2 dimer with a bond length of 2.86 A. 
Figure [3] shows both the symmetric and distorted config¬ 
urations, along with their band diagrams. The origin of 
the distortion is not due to the Jahn-Teller mechanism, 
because the symmetric state does not present degeneracy. 
As discussed above for the 2— state, the missing Te atom 
leaves four broken bonds with 0.5 electron each one. The 
energy is reduced when two Cd atoms loss 0.5 electrons 
and move to adopt a quasi-planar bonding structure with 
their three Te neighbors. The rules of covalent bonding 


are satisfied if these two Cd atoms distribute their 1.5 
electrons in three sp^ orbitals that make cr bonds with 
the Te sp^ orbitals. The other two Cd atoms, which also 
have 0.5 x 2 = 1 electron in dangling bonds, trap the 
electron released by the first Cd pair, getting closer to 
form a doubly-occupied covalent bond. This bond corre¬ 
sponds to the VB top in Fig. [3Kb), where the correspond¬ 
ing charge density (square of the wave function) is shown 
in Fig. SKa). 



FIG. 4. (Color online) Charge density isosurfaces of the defect 
state in the the distorted Te vacancy in the 0 (a) and 2— 
(b) charge states. Pink and cyan balls represent Te and Cd, 
respectively. The numbers indicate the Cd-Cd distances. 


For the 2— charge state, the total energy is apparently 
minimized by a D 2 d distortion leading to the formation 
of two perpendicular Cd dimers, as illustrated in Fig. |S] 
This case is instructive about the interplay between the 
different size effects. The band diagrams of Fig. [S] sug¬ 
gests that the band energy is minimized by a downshift of 
the defect band, that becomes occupied, while the disper¬ 
sive host-like conduction band gets empty. This mecha¬ 
nism lowers the total energy by 0.6 eV. However, the self- 
consistent energy must be corrected by subtracting the 
band filling energy, which is -0.87 eV for the undistorted 
configuration, but it is negligible for the distorted one. 
Therefore, the symmetric configuration has lower energy. 
In order to test this assertion, we have inserted the dis¬ 
torted configuration into a larger supercell composed of 
4x4x4 conventional unit cells, i.e. the SC512 model, 
which allows to use only the F point. With this setup, 
the band filling error is null. As expected, the distorted 
configuration in the SC512 relaxes to the non-distorted 
structure. Therefore, the distortion around the vacancy 
for the 2— state is a supercell-size effect driven by the 
conduction band filling in the SC64 and SC216 models. 
According to the above discussion, the distorted configu¬ 
ration is unstable within the^PBE approximation in the 
limit of dilute impurity. However, the quasiparticle cor¬ 
rections (see AppendixjB]) increase the total energy of the 
symmetric configuration, leaving the distorted configura¬ 
tion stable. The symmetric state has a PHS recognizable 
as the dispersive CB in Fig. [S] (a), which is occupiedby 
two electrons. On the other hand, in the distorted config- 
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electrons in the CB. Therefore the distorted configuration 
is stable. The calculation with the model SC512 con¬ 
firms this result. Moreover, the absence of conduction 
electrons, valence holes, and net charge, renders quasi¬ 
particle corrections unnecessary. 

The H- charge state has one electron occupying the 
CB in a PHS. The HSE06 calculation confirms that the 
PHS has lower energy than the localized defect state. 
It seems contradictory that for the 1-1- state, with one 
missing electron, the highest occupied level jumps into 
the CB, while for the 2+ state the highest occupied level 
returns to the VB. This behavior is due to the strong 
relaxation of the four neighbor Cd atoms into planar 
configurations CdTes, favoring a sp^-like hybridization. 
Within a sp^ shell structure, each Cd atom shares its 2 
valence electrons with 3 x 1.5 electrons coming from its 
three neighboring Te, resulting in 0.5 electron in excess 
of the sp^ shell. Adding the four Cd atoms, this sums 2 
electrons. The net effect is that two electronic levels are 
removed (holes). Therefore, the closed-shell structure is 
recovered in the charge state 2+, while for the state 1-1- 
there is one electron in the CB. Accordingly, for the 1-1- 
charge state, the formation energy must be corrected due 
to errors in both the conduction band filling and the VB 
and CB edges. For the latter, Ze = I and Zh =2 must be 
considered in Eq. IBll respectively. 


FIG. 5. (Color online) Stable geometries and band diagrams 
for Te vacancy in the 2— charge state withont (a), and with 
(b) a D 2 d distortion, nsing the SC216 snpercell. The solid 
black (green) lines are the bands of the defective (perfect) 
supercell, and the dashed lines represent the highest occupied 
level. Pink and cyan balls represent Te and Cd, respectively. 


uration, the two extra electrons occupy a localized state, 
shown in Fig. H^b), resulting in the fiat band of Fig. [5] 
(b). To assess the PBE gap error, we have computed the 
energy levels at the P point, using the HSE06 functional. 
In this calculation, only the P point was used to obtain 
the self-consistent charge density. This approximation is 
supported by comparing the corresponding PBE energies 
and wave functions using different k-point samplings, the 
P point and the 2x2x2 mesh. The effect of changing the 
PBE functional by the HSE06 functional is the shift of 
the PHS energy by 0.4 eV. In contrast, the defect level is 
shifted by only -0.07 eV. Additional GW calculations for 
the primitive cell, show that the CBM is 0.481 eV over the 
PBE value (see Appendix IB]) . Considering that the PHS 
is occupied by two electrons, the total energy of the sym¬ 
metric configuration must be corrected by 2 x 0.481 eV. 
However, for the distorted configuration the quasiparti¬ 
cle correction is quite small and can be neglected. As a 
consequence, the distorted configuration has the lowest 
energy. 

For the neutral state, there is no need of band filling 
corrections because there is neither holes in the VB nor 


C. Formation energy 

Figure [5] shows the formation energies of the Vxe de¬ 
fect performed with the SC64 supercell in several charge 
states. These energies do not contain the band filling 
and quasiparticle corrections. However, these corrections 
partially cancel each other, providing reasonable results. 
Thus, the differences with previous studiesSii^id^ can be 
easily understood. In Fig. [51 the negatively charged 
states appear stable for high values of the Fermi level 
because the energies of the excess electrons in the CB 
is heavily underestimated and, as discussed above, the 
quasiparticle correction (Eq. IbD) could not be applied. 
Recently, Biswas and Du^ using the hybrid PBEO func¬ 
tional tuned to fit the CdTe band gap, found no stable 
negative state. Well before, Wei and Zhangi^ using LDA, 
applied several approaches that led to a large band gap 
and no negative stable states. Later on, Du et al^ us¬ 
ing the LDA with quasiparticle correction, found that 
the stable states are 2-|-,0, and 2—. Their results were 
attributed to C 2 v and D 2 d distortions that reduce the 
energies of the states 0 and 2—, respectively. The forma¬ 
tion energies of Fig. [5] correspond to undistorted defects, 
obtained by relaxation of the structure after removal of a 
Te atom. Indeed, we found distorted configurations that 
lower the energy of the 0 and 2— states. According to the 
discussion of Section illl A[ we will analyze the distorted 
structures only for the larger SC216 supercell. 

Recent calculations using the HSE06 methodiSiii sug¬ 
gest that Vxe is stable only in the 2-\- charge state. These 






















6 


authors attribute the high stability of the charged state 
to the large Cd displacements around the vacancy. We 
think these authors may have missed the energy down¬ 
shift of the neutral and 2— states due to distortion. In 
any case, we consider unreliable the results obtained with 
the SC64 model. 
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FIG. 6. (Color online) Formation energies (AH) as a function 
of the Fermi energy (Ep) computed with a 64-atom supercell 
(SC64), not considering distortions. 


The formation energies, computed with the SC216 
model, are given as 
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Figure [ 7 ] shows the formation energy plots as function 
of the Fermi level. The quasiparticle corrections has been 
applied according to the Appendix [0 shifting the VBM 
by -0.470 eV and the CBM by 4-0.481 eV. Remarkably, 
the formation energy of the stable charge states 0 and 
2-|- are insensitive to the quasiparticle corrections. These 
corrections are not required for the neutral state, while 
for the 2-1- state the terms m and (IB2I) cancel each 
other. The net quasiparticle correction for the states 2—, 
1 —, and H- are 0.939, 0.951 and 0.951 eV, respectively. 
As this method gives a band gap of 1.55 eV, this is the 
maximum value for the Fermi level in Fig. [71 The ther¬ 
modynamic transition level e(2 -|- /O) = 0.69 eV. It is 
worth to note that almost the same value for e(2 -|- /O) is 
obtained with the model SC64 without band-filling and 
gap corrections (see Fig. 7). For the other charge states, 
the observed good agreement between Figs. [6] and [ 7 ] is 
due to error compensations between the gap error and 
the band-filling error. 

Some calculations have predicted negative formation 
energies for Vxe, mostly in Te-poor conditions but 
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FIG. 7. (Golor online) Formation energies as a function of the 
Fermi energy for the Te vacancy computed with a 216-atom 
supercell (SC216). 


also in Te-rich conditions^ Negative formation energies 
indicate a trend to high defect concentration and agglom¬ 
erate. High defect concentration need to be studied by a 
different methodology that correct this apparently wrong 
result. However, high concentration of Vxe in CdTe is 
not observed)^ (Ch. 5) It is worth to note that our cal¬ 
culations always predict positive formation energies for 
Vxe- 


IV. CONCLUSIONS 

In summary, the use of the 64-atoms supercell, as the 
basis to calculate defects in CdTe and other semiconduc¬ 
tors, is revisited. For the Te vacancy, Vxe, the band 
structure of some charge states presents unphysical fea¬ 
tures that are easily corrected by using larger supercells. 
These unphysical features, notably a dramatic alteration 
of the conduction band, makes uncertain to apply sim¬ 
ple quasiparticle corrections to the band edges. The for¬ 
mation energies and thermodynamic transition levels are 
slightly modified by these errors due to the partial com¬ 
pensation of the band-filling error and the gap underes¬ 
timation of GGA functionals. 

The band-filling error associated to the k-point sam¬ 
pling can drive the system to distorted configurations, as 
the case of . The calculations based on the 512-atoms 
cubic supercell SC512, with F—point sampling, are free 
of this error. However, the symmetric configuration turns 
out to have higher energy after quasiparticle (gap) cor¬ 
rections. In this sense, the band-filling error mimics the 
quasiparticle corrections, allowing to relax the system to 
the distorted configurations. For calculations when the 
SC512 model is not affordable, we suggest a careful anal¬ 
ysis of the band-filling errors in both distorted and undis¬ 
torted configurations. High level calculations like hybrid 
functionals, GW, or quantum Monte Garlo, are probably 
















7 


not affordable for such large supercells, and may need to 
be done for smaller supercells. In these cases, further op¬ 
timization of the distorted structures inserted in a large 
supercell with a single k-point and a low-cost GGA func¬ 
tional, are recommended. 

According to the literature, the Te vacancy seems to 
be the dominating donor and has strong influence on 
the electric and thermodynamical properties of CdTei^^— 
Our study shows that predictions based on calculations 
with the SC64 model are questionable, and additional 
efforts must be made, considering larger supercells like 
SG216 or SG512 in order to obtain reliable results. 
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Appendix A: Size corrections 

An isolated defect either populates the conduction 
band minimum (CBM) or create holes at the valence 
band maximum (VBM). However, in a finite supercell, 
when a k-point mesh is used, the electron or hole created 
by the point defect fills half a band, considering the spin 
degeneracy. The single-particle contribution to the total 
energy is affected by the band dispersion. The correction 
is to replace the energies of the defect band, either of 
electrons or holes, by the corresponding band extremum. 
For donors, the correction can be expressed as2i 

0(en,fe S.min)WkTln,k{s-n,k Gmin) j (Al) 

n,k 

where indexes n and k denote the bands and k-points, Wk 
are the k-point weights, en,k are the single particle ener¬ 
gies with minimum emim iln,k are the occupation num¬ 
bers, and 0 is the Heaviside function. The correction for 
acceptors is 

)wk{S - rin,k){e 

max 

n,k 

(A2) 

where S' = 1 (^ = 2) stands for spin polarized (unpolar¬ 
ized) calculations. Cmax is the maximum of the acceptor 
band. 

In VASP calculations, the average of the electrostatic 
potential is set to zero. Hence, the potential energy ref¬ 
erence is undetermined by an additive constant, which 


is different in supercells containing different number of 
atoms or different charge. The change has two sources: 
a) the presence of excess charge, b) the missing or appear¬ 
ance of extra atoms. For charged systems, the additive 
constant enters in the Hartree term of the total energy 
per cell: 


Eh = j e[no{r) - N{r) + An{r)]{v{r) + Av)d^r, 


where e is the electron charge, no(r) is the electron den¬ 
sity for a neutral system, An(r) is the excess electron 
density. iV(r) is the nuclear charge, and Av is the ad¬ 
ditive potential implicitly added to make zero the aver¬ 
age potential, with respect to the potential of the system 
without defect and neutral. The contribution of the ad¬ 
ditive potential to the Hartree energy is equal to 


AEh = J eAn{r)Avd^r 

An(r)d^r {eAv) 


{-q)AV. 


Here, q is the number of electrons donated 
f An(rd^r)/(—e), and AV = eAv is the change of 
electron potential energy. AV can be estimated by 
the core level shifts of atoms far from the defect with 
respect to the core levels in the calculation with the 
perfect crystal. The correction to the total energy due 
to potential alignment (p.a.) is 

AEp,a. = qAV. (A3) 

Periodic boundary conditions induce an spurious in¬ 
teractions of the carved defect with the periodic replicas. 
The leading terms can be corrected by the Makov-Payne 
formula 


AEm-p 


ctAiq^ 

2eVi/3 


W- + o(v-^'% 


where om is a Madelung constant, e is the dielectric con¬ 
stant, V is the supercell volume, and Q is a quadrupole 
moment. We have used a simplified version due to Lany 
and Zunger— 


AEl-z 


2 

3 2eVi/3' 


(A4) 


We have used the static dielectric constant e = 12.272, 
obtained from density functional perturbation theory cal¬ 
culations and including local field effectSi ^^i^^ 

In summary, the supercell size corrections are given by 
the sum of band-filling, potential alignment, and Lany- 
Zunger corrections 


AEgize — AEt.f. + AEpa. + AEi^-z- (A5) 







Appendix B: Quasiparticle corrections 


ence of the Fermi level 


Semilocal exchange-correlation functionals (LDA, 
GGA) cause a severe underestimation of the valence to 
conduction band gap, and this affects the formation en¬ 
ergies. Let us consider the gap error 


(B2) 

The total quasiparticle correction is 

AE,.p. = AE^^Hq) + AE^^\q). (B3) 


AEg = A-P - = AEc - AEv. 

The above equation established a partition of the gap er¬ 
ror in shifts applied to the valence and conduction band. 

For shallow acceptors, where active holes occupy per¬ 
turbed host states, the hole levels undergo the same shift 
as the valence band edge. This has been verified in 
quasiparticle calculationsFor shallow donors, the ac¬ 
tive electron levels are perturbed states of the conduction 
band, undergoing the shift. Including fractional occupa¬ 
tions, the shifts of single particle levels and total energies 
can be expressed as 


AF;«(g) = | 


Ze{X,q)AEc 

Zh{X,q){-AEv) 


for donors, 
for acceptors, 


(Bl) 


where Zg (zh) is the number of electrons introduced (re¬ 
moved) by the defect near the CBM (VBM). 

A second quasiparticle correction modifies the refer- 


With our computational setup with the PBE func¬ 
tional, the VBM and CBM energies are 1.789 and 2.388 
eV, respectively. For quasiparticle corrections we have 
considered the generalized DFT with the hybrid func¬ 
tional HSE06i^ and GW calculations^ including the 
spin-orbit coupling, for the CdTe primitive cell, using 
a 8 X 8 X 8 k-point grid. The HSE06 VBM and CBM 
energies are 1.773 and 2.977 eV, respectively. With this 
approximation the band gap is 1.204 eV, which is still 
lower than the experimental value. The single-shot per¬ 
turbative GqWo calculatioii^i based upon the HSE06 cal¬ 
culation provides a band gap of 1.55 eV, and downshift 
of the VBM by 0.454 eV below the HSE06 value. The 
latter HSE06 and GqWo were performed with PAW po¬ 
tentials adapted to GW calculations. The used cutoff for 
the response function was 169.4 eV. The integrals defin¬ 
ing the self-energies were evaluated using 96 frequencies 
between 0 and 288 eV /h. The number of (occupied plus 
empty) bands was carried up to 1600, and the VBM was 
obtained by extrapolation to infinite bands. 
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